Ref. SISSA 26/2006/EP 
May 2006 
hep-ph/0605151 



Charged Lepton Decays li — ► lj + 7, Leptogenesis CP- Violating 
Parameters and Majorana Phases 

S. T. Petcov 1 and T. Shindou 

Scuola Internazionale Superiore di Studi Avanzati, and 
Istituto Nazionale di Fisica Nucleare, 1-3401 4 Trieste, Italy 



Abstract 

We analyse the dependence of the rates of the LFV charged lepton decays \x — > e + 7, 
r — > e+7, r — ► /U+7 (Zj — > /j+7) and their ratios, predicted in the class of SUSY theories 
with see-saw mechanism of i^-mass generation and soft SUSY breaking with universal 
boundary conditions at the GUT scale, on the Majorana CP-violation phases in the 
PMNS neutrino mixing matrix and the "leptogenesis" CP-violating (CPV) parameters. 
The case of quasi-degenerate in mass heavy Majorana neutrinos is considered. The 
analysis is performed for normal hierarchical (NH), inverted hierarchical (IH) and quasi- 
degenerate (QD) light neutrino mass spectra. We show, in particular, that for NH and 
IH z^-mass spectrum and negligible lightest neutrino mass, all three li — ► lj + 7 decay 
branching ratios, BR(li — * lj +7), depend on one Majorana phase, one leptogenesis 
CPV parameter and on the 3-neutrino oscillation parameters; if the CHOOZ mixing 
angle #13 is sufficiently large, they depend on the Dirac CPV phase in the PMNS matrix. 
The "double ratios" R(21/31) ~ BR{u -> e + j)/BR(t -> e + 7) and R(21/32) ~ 
BR(fi — > e + 7) / BR(t — > fi + 7) are determined by these parameters. The same 
Majorana phase enters into the NH and IH expressions for the effective Majorana 
mass in neutrinoless double beta decay, <m>. 



1 Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, 
Bulgaria. 
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1 Introduction 



The experiments with solar, atmospheric, reactor and accelerator neutrinos [1-5] have 
provided during the last several years compelling evidence for the existence of non-trivial 
3-neutrino mixing in the weak charged-lepton current (see, e.g., [6]): 

3 

vvl = ^ Ul i ^' L ' 1 = e ' ^> r ' W 

where vil are the flavour neutrino fields, i/^x is the field of neutrino Vj having a mass rrij 
and U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [7], U = ?7pmns- 
The existing data, including the data from the 3 H /3-decay experiments [8] imply that the 
massive neutrinos Vj are significantly lighter than the charged leptons and quarks: rrij < 2.3 
eV (95% C.L.) '. 

The existence of the flavour neutrino mixing, eq. implies that the individual lepton 
charges, L/, I = e,\i,r, are not conserved (see, e.g., [11]), and processes like \x~ — > e~ + 7, 
/i~ — > e~ + e + + e _ , t~ — ► e~ + 7, r~ — > /i _ + 7, /i~ + (A, Z) — > e~ + (A, Z) , etc. should take 
place. Stringent experimental upper limits on the branching ratios and relative cross-sections 
of the indicated |AL/| = 1 decays and reactions have been obtained [12-14] (90% C.L.): 

BR(/i -> e + 7) < 1.2 x KT 11 , BR(/x -> 3e) < 1.2 x 10" 12 , 
BR(r -> /i + 7) < 6.8 x 10~ 8 , R(/i" + Ti -> e~ + Ti) < 4.3 x 10~ 12 . <_) 

Future experiments with increased sensitivity can reduce the current bounds on BR(/i — > 
e + 7), BR(r — > /x + 7) and on R(/i~ + (A, Z) — >• e _ + (A, Z)) by a few orders of magnitude 
(see, e.g., [15]). In the experiment MEG under preparation at PSI [16] it is planned to reach 
a sensitivity to 

BR(ii^e + 7) ~ (1(T 13 - 1(T 14 ) . (3) 

In the minimal extension of the Standard Theory with massive neutrinos and neutrino 
mixing, the rates and cross sections of the LFV processes are suppressed by the factor [17] 
(see also [18]) (m^/M^y) 4 < 6.7 x 1CT 43 , M w being the W ± mass, which renders them 
unobservable. It was shown in [19] that in SUSY theories with see-saw mechanism of 
neutrino mass generation 2 [21] and soft SUSY breaking with universal boundary conditions 
at a scale Mx above the right-handed (RH) Majorana neutrino mass scale Mr, Mx > Mr 
the rates and cross sections of the LFV processes can be strongly enhanced and can be within 
the sensitivity of presently operating and future planned experiments (see also, e.g., [22-31]). 
As is well-known, the see-saw mechanism of neutrino mass generation [21], provides a very 
attractive explanation of the smallness of the neutrino masses and - through the leptogenesis 
theory [32], of the observed baryon asymmetry of the Universe. 

One of the basic ingredients of the see-saw mechanism is the matrix of neutrino Yukawa 
couplings, Y u . Leptogenesis depends on Y u as well. In the large class of SUSY models 

^^More stringent upper limit on rrij follows from the constraints on the sum of neutrino masses obtained 
from cosmological/astrophysical observations, namely, the CMB data of the WMAP experiment combined 
with data from large scale structure surveys (2dFGRS, SDSS) [9]: £V rrij < (0.7-2.0) eV (95% C.L.), where 
we have included a conservative estimate of the uncertainty in the upper limit (see, e.g., [10]). 

2 An integral part of the see-saw mechanism are the right-handed (heavy) Majorana neutrinos [20]. 
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with see-saw mechanism and SUSY breaking mediated by flavour-universal soft terms at a 
scale M x > M R we will consider, the probabilities of LFV processes also depend strongly on 
Y u (see, e.g., [23,24]). The matrix Y u can be expressed in terms of the light neutrino and 
heavy RH neutrino masses, the neutrino mixing matrix £/pmns> an d an orthogonal matrix 
R [23]. Obviously, Y u depends on the Majorana CP- violation (CPV) phases in the PMNS 
matrix £/pmns [33]. In the case of negligible flavour effects in leptogenesis [34], successful 
leptogenesis is possible only if R is complex (see, e.g., [35]). For M R < 10 12 GeV, the 
flavour effects in leptogenesis can be substantial [34] (see also [36,37]). Due to the latter, 
leptogenesis can take place even for real R [34,38], but only if R ^ 1. In this case the Dirac 
and/or the Majorana CPV phases in the PMNS matrix play the role of the CPV parameters 
responsible for the generation of the baryon asymmetry of the Universe. Thus, in this case 
there is a direct link between the low-energy leptonic CP- violation and the generation of the 
baryon asymmetry of the Universe [38]. 

It was shown in [26] that if both the light and the heavy Majorana neutrino mass spectra 
are quasi-degenerate (QD), the rates of LFV decays \i — > e + 7, r — > e + 7 and r — > 
/i + 7, predicted in the class of SUSY theories of interest, can be strongly enhanced by the 
leptogenesis CP-violating (CPV) parameters in the complex matrix R, with respect to the 
rates predicted for real R ^ 1 or for R = 1. The indicated LFV decay rates were also 
noticed in [26] to depend for complex 1 on the Majorana CPV phases in L/pmns- This 
dependence was investigated recently in [30] by taking into account the effects of the phases 
in the renormalisation group (RG) running of the light //-masses rrij and of the mixing angles 
in L/pMNS- It was found [30] that the Majorana phases can affect significantly the predictions 
for the fj, — > e + 7 and r — * e + 7 decay rates. 

In the present article we extend the analyses performed in [26,30] to the cases of normal 
hierarchical and inverted hierarchical light neutrino mass spectra. We investigate also in 
greater detail the case of QD spectrum. More specifically, working in the framework of 
the class of SUSY theories with see-saw mechanism and soft SUSY breaking with flavour- 
universal boundary conditions at a scale Mx > Mr, we study in detail the dependence of the 
rates of charged lepton flavour violating (LFV) radiative decays /i ^ e + 7, r ^ e + 7 and 
t — > /i + 7, on the Majorana CPV phases in {/pmns an d on the leptogenesis CPV parameters 
in the complex orthogonal matrix R. The case of quasi-degenerate (QD) in mass heavy 
RH Majorana neutrinos is considered. It is well-known that in the case of heavy Majorana 
neutrinos with QD mass spectrum (i.e., negligible splitting between the masses), the rates 
of LFV radiative decays of interest do not depend on the matrix R 7^ 1 if R is a real matrix 
(see, e.g., [26]). Our analysis is performed under the condition of negligible RG effects for 
the light neutrino masses rrij and the mixing angles and CP- violation phases in t/pMNS- The 
RG effects in question (see, e.g., [30,39] and the references quoted therein) are negligibly 
small in the class of SUSY theories we are considering in the case of hierarchical (normal 
or inverted) Vj mass spectrum. The same is valid for QD Vj mass spectrum provided the 
SUSY parameter tan j3 is relatively small, tan j3 < 10, tan (3 being the ratio of the vacuum 
expectation values of the up- and down- type Higgs doublet fields in SUSY extensions of the 
Standard Theory. For the three types of light neutrino mass spectrum, we investigate also 
the predictions for the ratios of the rates of fj, — > e + 7 and r — > e + 7, and of fi — * e + 7 
and r — > /i + 7, decays. In a large region of the relevant SUSY parameter space these 
two ratios are independent of the SUSY parameters and are determined completely by the 
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neutrino mixing angles, Majorana and Dirac CPV phases, leptogenesis CPV parameter(s) 
and, depending on the type of the neutrino mass spectrum - hierarchical or quasi-degenerate, 
by the neutrino mass squared differences Am^ and Am^ or the absolute neutrino mass. A 
study of the predictions for the two LFV decay rate ratios was performed recently in ref. [40] . 
In [40], however, only the case of zero leptogenesis CPV parameters (i.e., of real matrix R) 
and of zero Majorana CPV phases in {/pmns was investigated. 



2 Neutrino Mixing Parameters from Neutrino Oscillation Data 

We will use the standard parametrisation of the PMNS matrix ?7 PM ns (see, e.g., [41]): 



PMNS 



/ C12C13 S12C13 s l3 e lS \ 

-S12C23 - Ci 2 S 2 3Sl3e l<5 C12C23 - s 12 s 2 3Si3e t5 S23C13 
V S12S23 — c 12 c 23 s 13 el5 ~~ C12S23 — s 12C23 s 13 eJ<5 C23 c 13 / 



diag(l,e i 2,e*' 3 2 f ) , 



where c. 



(4) 

cos%, Sij = sm9ij, the angles % = [0,7r/2], 5 = [0,2n] is the Dirac CP- 
violating phase and a and {3m are two Majorana CP-violation phases [33,42]. One can 
identify the neutrino mass squared difference responsible for solar neutrino oscillations, Am 2 & , 



with Ar/ijj = m\ 
the dominant v„ 



m\, Am & 



Am 2 21 



\Am\\ 



\ Am li\ 



v T [i/ f 
\Am 2 32 \ 



> 0. The neutrino mass squared difference driving 
— > u T ) oscillations of atmospheric (P M ) is then given by 
Am^. The corresponding solar and atmospheric neutrino 



mixing angles, & and 9\, coincide with 9 12 and #23, respectively. The angle #13 is limited 
by the data from the CHOOZ and Palo Verde experiments [43]. 

The existing neutrino oscillation data allow us to determine Am^, jAm^l, sin 2 # 12 and 
sin 2 2#23 with a relatively good precision and to obtain rather stringent limits on sin 2 # 13 
(see, e.g., [44,45]). The best fit values and the 95% C.L. allowed ranges of Am^, sin 2 On, 
(Am^l and sin 2 2^23 read 3 : 



Am 2 21 = 8.0 x 10~ 5 eV 2 , sin 2 21 = 0.31 , 
Am 2 x = (7.3 - 8.5) x lO" 5 eV 2 , sin 2 12 = (0.26 - 0.36) , 



(5) 



IATO31I 



(6) 



2.2xlO" 3 eV 2 , sin 2 2 23 = 1.0 , 
1.7 - 2.9) x 10" 3 eV 2 , sin 2 2 23 > 0.90. 
A combined 4 3-u oscillation analysis of the solar neutrino, KL and CHOOZ data gives [44] 

sin 2 fl 13 < 0.027 (0.044), at 95% (99.73%) C.L. (7) 

The neutrino oscillation parameters discussed above can (and very likely will) be measured 
with much higher accuracy in the future (see, e.g., [6]). 

3 The data imply, in particular, that maximal solar neutrino mixing is ruled out at ~ 6a; at 95% C.L. one 
finds cos 200 > 0.26 [44], which has important implications [46]. 

4 Using the recently announced (but still unpublished) data from the MINOS experiment [47] in the 
analysis leads to somewhat different best fit value and 95% allowed range of jAm^l [48]: jAm^l = 2.5 x 
1(T 3 eV 2 and \Am 2 31 \ = (2.2 - 2.9) x 1CT 3 eV 2 . 
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The sign of Am\ = Am\ x , as it is well known, cannot be determined from the present (SK 
atmospheric neutrino and K2K) data. The two possibilities, AmL 32 j > or Am^^ < 
correspond to two different types of z/-mass spectrum: 

- with normal ordering (hierarchy) mi < m 2 < m 3 , Am\ = Ara^ > 0, and 

- with inverted ordering (hierarchy) m 3 < mi < m 2 , Am\ = Am^ < 0. 

Depending on the sign of Am^, sgn(Am^), and the value of the lightest neutrino mass, 
min(mj), the v-m&ss spectrum can be 

- Normal Hierarchical: uii^ m 2 <^m 3 , m2=(Amg)^~ 0.009 eV, m3=|Am^|^~ 0.05 eV; 

- Inverted Hierarchical: m 3 <C mi < m 2 , with mi j2 = lAm^s ~ 0.05 eV; 

- Quasi-Degenerate (QD): mi = m 2 = m 3 = m, m 2 ^> |Am^|, m ^ 0.10 eV. 

The sign of Am\ x = Am\ 2 , which drives the dominant atmospheric neutrino oscillations, 
can be determined by studying oscillations of neutrinos and antineutrinos, say, — > v e and 
&n v e , i n which matter effects are sufficiently large. This can be done, e.g., in long-baseline 
//-oscillation experiments (see, e.g., [49]). Information about sgn(Am| 1 ) can be obtained also 
in atmospheric neutrino experiments by studying the oscillations of the atmospheric and 

which traverse the Earth [50]. 

As is well-known, the theories employing the see-saw mechanism of neutrino mass gen- 
eration [21] of interest for our discussion, predict the massive neutrinos Uj to be Majorana 
particles. Determining the nature of massive neutrinos is one of the most formidable and 
pressing problems in today's neutrino physics (see, e.g., [6,51]). If it is established that the 
massive neutrinos Uj are indeed Majorana fermions, getting information about the Majorana 
CP- violation phases in t/pMNS; would be a very difficult problem. The oscillations of flavour 
neutrinos, V\ —>■ vy and V\ — > vy , 1,1' = e,fi,r, are insensitive to the Majorana CP- violation 
phases a and (3m [33,52]. The only feasible experiments that at present have the potential 
of establishing the Majorana nature of light neutrinos Uj and of providing information on 
the Majorana CP-violation phases in {/pmns are the experiments searching for the neutri- 
noless double beta ((/5/3)o^-) decay, {A, Z) — > {A, Z + 2) + e~ + e~ (see, e.g., [11,51,53]). 
The {(3(3) Qv -decay effective Majorana mass, <m> (see, e.g., [11]), which contains all the 
dependence of the {(3(3) o-decay amplitude on the neutrino mixing parameters, is given by 
the following expressions for the normal hierarchical (NH), inverted hierarchical (IH) and 
quasi-degenerate (QD) neutrino mass spectra (see, e.g., [53]): 



\<m>\ 



Am 2 21 sin 2 6 u e ia + J Am 2 31 sin 2 9 13 e i0M 



, mi«m 2 «m 3 (NH), (8) 



|<m>| = \J Am\ 3 |cos 2 6»i 2 + e ia sin 2 6»i 2 | , m 3 < m x < m 2 (IH), (9) 

|<m>| ^ m |cos 2 #i 2 + e ia sin 2 fl 12 | , m i)2 ,3 = m > 0.10 eV (QD). (10) 

Obviously, |<m>| depends strongly on the Majorana CP- violation phase(s) 5 ; the CP- 
conserving values of {a — (3m) = 0, ±7r {a = 0, ±7r) [54], in particular, determine the range 

5 We assume that the fields of the Majorana neutrinos Vj satisfy the Majorana conditions: C(vi^) T = 
and C(y-i) T = e~ l2S vz, where C is the charge conjugation matrix. With the parametrisation we are employing 
for C/pMNS, eq. (@J), the effective Majorana mass |<m>| does not depend on the Dirac CP-violation phase 
8 as a consequence of the presence of the phase factor e~ l2S in the Majorana condition for the field ^3. 
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of possible values of |<m>| in the case of NH (IH, QD) spectrum. If the (/3/3)o I y-decay is 
observed, the measurement of the (f3/3) o u -dec&y half-life combined with information on the 
absolute scale of neutrino masses (or on min(mj)), might allow to significantly constrain the 
Majorana phase a [41,55,56], for instance. 

3 The See-Saw Mechanism, Neutrino Yukawa Couplings, and 
LFV Decays k — > lj + 7 

In the minimal supersymmetric standard model with RH neutrinos Nj and see-saw 
mechanism of neutrino mass generation (MSSMRN) we consider it is always possible to 
choose a basis in which both the matrix of charged lepton Yukawa couplings, Ye, and the 
Majorana mass matrix of the heavy RH neutrinos, Mn, are real and diagonal. We will work 
in that basis and will denote by Dn the corresponding diagonal RH neutrino mass matrix, 
Dn = diag(Mi, M 2 , M 3 ), with Mj > 0. We will consider in what follows the case of QD 
heavy Majorana neutrinos: Mi = M 2 = M 3 = Mr. It will be assumed that the splittings 
between the masses of the heavy Majorana neutrinos are sufficiently small, e.g., that they 
are of the order of those considered in [26] . The existence of sufficiently small (but nonzero) 
splittings between the masses of the heavy Majorana neutrinos Nj is indeed a necessary 
condition for the successful (resonant) leptogenesis to take place. The requisite small mass 
splittings can be generated, e.g., by renormalisation group effects [57]. However, the mass 
splittings under discussion, |Mj — Mj\ <C Mj, Mj, i ^ j — 1,2, 3, do not play any significant 
role in the predictions for the rates of the decays U — > lj + 7, which is the main subject of 
the present study. The heavy Majorana neutrino mass Mr will standardly be assumed to 
be smaller than the GUT scale M GUT ~ 2 x 10 16 GeV. 

In the class of theories of interest, the branching ratio of the /, — > lj + 7 decay has the 
following form (in the "mass insertion" and leading- log approximations, see, e.g., [22,24,28]): 

2 



BR(k -> Iji) 



3 



cm 



"total (0 G 2 F m\ 



s 



(3 + al)ml 



Y(Y^) ik ln^ Y kl 
^ { u > M k v 

k 



tan 2 /3, (11) 



where i ^ j = 1,2, 3, l\, l 2 , h = e, //, r, m and A = a m are the universal SUSY breaking 
scalar masses and trilinear scalar couplings at Mx > Mr, mj represents SUSY particle mass 
(see further), tan/3 is the ratio of the vacuum expectation values of up-type and down-type 
Higgs fields and Y u = Y u (Mr) is the matrix of neutrino Yukawa couplings evaluated at Mr. 
The matrix Y u can be parametrised as [23] 

Y„(Mr) = — a/d^ R a/D^ U 1 " = — a/Mr R v // D^ U 1 " . (12) 

Here v u = wsin/3, where v = 174 GeV, R is a complex orthogonal matrix 6 R T R = 1, 
= diag(mi, m2, W3), 77112,3 > being the light neutrino masses 7 and U is the PMNS 
matrix. 



6 Equation (|12|l represents the so-called "orthogonal" parametrisation of Y„. In certain cases it is more 
convenient to use the "bi- unitary" parametrisation [27] Y„ = U^ Y^ lag U^, where Ul,r are unitary matrices 
and Y^ lag is a real diagonal matrix. The orthogonal parametrisation is better adapted for our analysis and 
we will employ it in what follows. 

7 To be more precise, we can have mm(irij) = 0. 
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In what follows we will consider the case when the RG running of rrij and of the pa- 
rameters in t/pMNS from approximately Mz ~ 100 GeV, where they are measured, to M R 
is relatively small and can be neglected. This possibility is realised in the class of theories 
under discussion for sufficiently small values of tan/3 and/or of the lightest neutrino mass 
mm(rrij), e.g., for tan/? < 10 and/or mm(rrij) ^ 0.05 eV (see, e.g., [30,39]). Under the 
indicated condition, D„ and U in eq. (J12)) should be taken at the scale ~ Mz, at which the 
neutrino mixing parameters are measured. 

It was shown in [28] that in a large region of the relevant soft SUSY breaking parameter 
space, the expression 

m| ~ 0.5 ml m\i 2 (ml + 0.6 my 2 ) 2 , (13) 

mi/2 being the universal gaugino mass at Mx, gives an excellent approximation to the results 
obtained in a full renormalisation group analysis, i.e., without using the leading- log and the 
mass insertion approximations. For values of the soft SUSY breaking parameters implying 
SUSY particle masses in the range of few to several hundred GeV, say, m = mi/ 2 = 250 
GeV, A = a m = —100 GeV, we have: 

BR(k -> Ijj) ^ 9.1 x 10~ 10 KY+LY^I 2 tan 2 /3 , (14) 

where L ^ \n(M x /M R ). Since tan 2 (3 > 10, eq. JUJ) implies that if indeed the SUSY particle 
masses do not exceed several hundred GeV, the quantity |(YjLYy) 2 i| has to be relatively 
small. This is realised for, e.g., M R ^ 10 12 GeV. 

As follows from eqs. (fTTj) and (fTlj) and was widely discussed , in the case of soft SUSY 
breaking mediated by soft flavour-universal terms at Mx > Mr, the predicted rates of LFV 
processes such as /i — > e + 7 decay are very sensitive to the off-diagonal elements of 

Yl(M R )Y u (M R ) = \ Rt B N R v^Ut S ^ Rt R ^D v \]^ . (15) 

K K 

It is well-known that in the theories with see-saw mechanism, leptogenesis depends on 
Y V (M R ) and thus on R. In the case of negligible flavour effects [34], the dependence of 
interest is realised through the product [35] 

Y u (M R )Yl(M R ) = \ a/Dn R Rt = ^ R R f . (16) 

v v 

In this case successful leptogenesis can take place only if R 7^ 1 is complex. If M R ^ 
10 12 GeV, flavour effects in leptogenesis can be significant and leptogenesis can proceed 
successfully even for real 1 (see, e.g., [38]). It follows from eqs. (fTTj) and (fTKj) . however, 
that in the case of QD in mass heavy RH Majorana neutrinos of interest, the predicted rates 
of LFV decays fi — > e + 7, etc. are independent of the orthogonal matrix R if R is real. 

In what follows we will consider (R)* 7^ R and will use the parameterizations of R 
proposed in [26]: 

R = O e iA . (17) 
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Here O is a real orthogonal matrix and A is a real antisymmetric matrix, (A) T = —A, 



A 




i, b, c being real parameters. The following representation of e tA proves useful [26]: 



eiA = i _ coshr-l A2 + i Sinhr Ai 



with r = y/a 2 + b 2 + c 2 . The requirement of successful leptogenesis in the case of QD light 
and heavy Majorana neutrino mass spectra and negligible flavour effects [34] implies [26] 
that abc ^ and that none of the parameters |a|, \b\ and |c| can be exceedingly small: 
|a6c| ~ (10~ 6 — 10~ 4 ). One also finds from the condition that Yukawa couplings should have 
moduli which do not exceed ~ 1 that typically r < 1 [26]. 

The parametrisation given in eq. (|17|) is particularly convenient in the analysis of the case 
of QD heavy Majorana neutrinos. We will consider a range of values of the parameters a, b, c 
determined by 10~ 4 \a\, \b\, |c| 0.10. Equations (fTTj) and (fTo]) imply that for QD heavy 
Majorana neutrinos we can set O = 1 and use R = e tA in the calculation of BR(Zj — > lj + 7) 
without loss of generality. Results for BR(7j — > lj + 7) in the case of real R / 1 can be 
obtained by formally replacing iA by in the expressions for BR(/j — > lj + 7) derived using 
R = e iA . 

4 The LFV Decays U —> lj + 7 and Majorana Phases 

In the case under discussion Mi = M 2 = M 3 = Mr and the matrix of neutrino Yukawa 
couplings has the form ~Y U = w ^^ Oe tA ^'D L ,\J^ . The off-diagonal elements of Y^Y^ of 
interest do not depend on O and satisfy (Y^Y^)^- = (Y^Y^)^. To leading order in small 
quantities they are given by (see also [26,30]) 



M R 



(YjY u ) 12 = A 21 C 23 C 12 S U + A 3 iS23Sl 3 e l& 

a^m x m 2 {c 23 (c\ 2 e~ l * + s\ 2 e^) + 2i S13C12S12S23 e' t5 sin 
+ 6 x /m 1 m 3 S23 {c 12 e~ 1 ^ - Si 3 s l2 e l{L ^~' s) ) 
+ cy/m 2 m 3 i^s 23 s 12 e l ^~ - c 23 s 13 c 12 e~ t( ^^ +d) J + 0{s{ 3 



(20) 
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u)13 



Mr 



A21S23C12S12 + A 3 iC23Si 3 e 15 

f /O ■ a n ■ ol v ■ r Ot \ 

a y /m 1 m 2 y-s 23 (c 12 e 2 + Sl2 e 2 ) + 2i s 13 c 12 s 12 c 23 e sin -J 



+ 6 v /rnlm 3 " (ci 2 c 23 e~ 



1 2 



Sl3Sl2S 23 e 2 



+ c v /m 2 m 3 si 2 c 23 e 2 + si 3 Ci 2 s 23 e 2 ^ ■- 



+ 0(4) 



+ 0(r\ s{ 3 ) , 



(21) 



(^"J^)23 — A 31 S2 3 C2 3 



Mr 



^m 1 m 2 (-2ic 12 s 12 c 23 s 23 sin ^ + s i3 [c 2 12 (s 2 2Z e 



+W23 e 



2 „i(f+«) 



c 23 e 



)]) 



+ b y /m 1 m 3 -s 12 {s 23 e 2 + c 23 e 2 ) + si 3 ci 2 c 23 s 23 2zsm(— d) 



/ / 2 ,■ "—Pm > 

+ Cy/m 2 m 3 ci 2 (c 23 e 2 + s 23 e 



S13S12C23S23 2isin 



13m 
2 



+ 0(s 



13) 



where 



M h 



M, 



\rrii 



mi 



R 



Aw| 
rrii + rrij 



(22) 



(23) 



Equations ()2U )) -([22 )) are valid for any of the possible types of light neutrino mass spectrum. 
The corresponding expressions for real 1 can be obtained by formally setting the three 
leptogenesis CPV parameters a, b and c to in eqs. (j2THl - (J22Jl . 

In what follows we will concentrate on the case of complex R / 1 and will call the 
real quantities a,b,c "leptogenesis CP-violation (CPV) parameters". The results in eqs. 
(J20J) -()22 J) imply that in the absence of significant RG effects, the "double" ratios 



R(21/31) 



BR(/i -> e + 7) 
BR(r -»• e + 7) 



BR(r 



ev T v P 



R(21/32) 



BR(/i -> e + 7) 
BR(r /i + 7) 



BR(r 



ev T v e 



(24) 



depend in the region of validity of eqs. (|TT|) and (fT3*|) in the relevant SUSY parameter space, 
on the neutrino masses rrij, mixing angles 8u, 8 23 , 9\ 3 and Majorana and Dirac CP- violation 
phases a, [3m and 5 at ~ M z , as well as on the leptogenesis CP-violating (CPV) parameters 
a, b and c. The dependence of the Dirac phase 5 can be significant only if the CHOOZ 
angle 9\ 3 is sufficiently large. The case of real matrix R ^ 1 corresponds to a = 0, b = 
and c = 0. Although the general expressions for (YjY u )ij, i ^ j, eqs. (j2*U|l - (j2*2*Jl . include 
several terms, there are few physically interesting cases in which the expressions simplify 
and the dependence on the Majorana CP- violation phase(s) and/or on the leptogenesis CPV 
parameters is prominent. 
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4.1 Normal Hierarchical Neutrino Mass Spectrum 

If the neutrino mass spectrum is of the normal hierarchical (NH) type and m% is negligi- 
bly small, i.e., \a\ y /m 1 m 2 , \b\^/mim 3 -C \c\yjm 2 m 3 , the quantities (YjY v )ij, i ^ j, depend, in 
particular, on the same Majorana phase difference (a — (3m) on which the effective Majorana 
mass, eq. (jSJ), depends, on the Dirac CPV phase 5 and one leptogenesis CPV parameter, 
c. The terms oc c^fm 2 rn 3 s V3 e' %& give always subdominant contributions in | (Y$Y V ) 12,13 1- 
For S13 <C tan ^12 ~ 0.65 they are negligible. In this case the expressions for |(yj5^)i2,i3| 
simplify: 

\(YXm = |c 23 P NH + 323 Q NH | 2 , (25) 
\(YXm = I--2S P NU + c 23 Q NH | 2 , (26) 



i/ 1 

u 



where 

P NH = (Am^ Cl2Sl2 , (27) 
Q NH = (A m ^)^ 13 e- tf + i 2c (Am^Am^i s 12 e 1 ^ . (28) 

The double ratio R(21/31) is determined completely by the solar and atmospheric neutrino 
oscillation parameters Am 21 , #12 and Am 31 and #23; by the CHOOZ angle 8 13 and by the 
Majorana and Dirac CPV phases (a — /3m) and 5 and by the leptogenesis CPV parameter c: 

|c 23 P NH + S 23 Q NR \ 2 

R(21/31) = V — L . (29) 

I-S23 P NH + C 23 Q NH | 2 1 ^ 

where P NH and Q NH are given by eqs. (|2*7j) and (|2*5|) . It follows from eqs. (j2ZJ), (|28)i and 

that if (a — (3m) = and 5 = ±vr/2, ±37r/2, we would have 

2 IpNHI 2 , 2 InNHI 2 

Rf21/311 = 23 ' ' 23 l_£ I /'oqx 

s 2 3 |pNH|2 + c 2 3 |gNH|2 " 

For the best fit value sin 2 29 23 = 1 we get R(21/31) = 1 independently of the value of the 
leptogenesis CPV parameter c, although the corresponding branching ratios BR(|U — > e + 7) 
and BR(r —>■ e + 7) can exhibit strong dependence on c. If #23 differs somewhat from 7r/4, 
the dependence of R(21/31) on c will, in general, be relatively mild. For |p NH | 2 ^> |Q NH | 2 
(|P H | 2 <C |Q NH | 2 ), however, the dependence of R(21/31) on c will be negligible even if 
# 23 ^ tt/4 and we would have R(21/31) = cot 2 6 23 (tan 2 6 23 ). 

For |c| ^5 0.1 and s 13 having a value close to the existing (3<r) upper limit s 13 = 0.2, 
the term oc c in Q NR , eq. (|28|). gives practically negligible contributions in | (Y^Y^) 12,13 1 , 
while P NH , eq. (J27|) . gives a subleading but non-negligible contribution. Correspondingly, 
the double ratio R(21/31) exhibits in this case a significant dependence on the Dirac phase S 
and is essentially independent of the leptogenesis CPV parameter c and the Majorana phase 
(a-(3 M ) (Fig. 1). 

If S13 ^ 0.1, but S13 is not much smaller than a/ Am 21 ci 2 si 2 / a/ Am^ ~ 0.07, the branch- 
ing ratios BR(/i — > e + 7) and BR(r — > e + 7) still depend on the CHOOZ mixing angle #13 
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and the phase 5. For 0.01 ^ |c| ^ 0.10, BR(/i — > e + 7) and BR(r — > e + 7) can exhibit 
significant dependence also on the Majorana phase (a — (5m)- The dependence of the double 
ratio R(21/31) on (a — (3 m) and 8 can be very strong due to possible mutual compensation 
between P NH and Q NH (see eq. (JUJ)). For {a-(3 M ) = 0, s i3 = 0.10 and sufficiently small |c|, 
for instance, we can have R(21/31) ~ 10~ 2 or R(21/31) ~ 10 2 depending on whether 5 = 71 
or 5 = 0; for (a - (3 M ) = vr and |c| ~ 0.03, R(21/31) can have a value R(21/31) ~ 10~ 3 or 
R(21/31) ~ 10 3 , respectively (Fig. 1). 

For rather small values of s 13 , namely, S13 <C a/ Am 2 1 ci 2 s 12 / ^/Am^, the dependence of 
1(1^1^)12,13! 011 ^is 6 "* 5 is insignificant and can be neglected. Under the latter condition we 
also have a/ Am^^ <C a/ Am^s^. The effective Majorana mass in (/3/3)o I /-decay is given 
correspondingly by |<m>| NH = a/ Am^s^. The quantities P NH and Q NH can be written 
as: 

P NH S (|<m>| NH )5(Am 2 1 )ici 2 , Q NH = 2 2c (|<m>| NH )l(Am 2 1 )3 . (31) 

Thus, in this case BR(/i — * e + 7) oc |<m>| NH and BR(r — ► e + 7) oc |<m>| NH . Given 
Am^i, Am 2 i, 6*12 and 623, the ratio R(21/31) is determined by the Majorana phase difference 
(a — (3m) and the leptogenesis CPV parameter c: 



R(21/31) 



(Am2i)4Ci2 cot 9 2 3 + i 2c (Am^)* e 



. <X-B M 

1 2 



2 



(Am2i)4Ci2 — i 2c (Amy! e 4 2 A/ cot 623 



(32) 



Obviously, for (a — (3m) — 0, we have R(21/31) = 1. If, however, (a — (3m) — ±tt, the double 
ratio R(21/31) can depend strongly on the value of |c|, provided |c| ^ 0.05. For 2|c| ~ 
(Am^i/Am^)* = 0.42, the two terms in the numerator (denominator) of the expression for 
R(21/31) can compensate (partially) each other and one can have R(21/31) ~ (10~ 3 — 10~ 2 ) 
or R(21/31) ~ (10 3 — 10 2 ) depending on the sign of c (Fig. 1). 

If |c| is relatively small, 2\c\ < (Am^/Am^)* = 0.42, BR(/i -> e+7) and BR(r -> e+7) 
are practically independent of c and (a — (3m)- This case was analised recently in [30]. If, 
for instance, Si 3 <C a/ Am^c^s^/ \/ Am^, we find from eq. (|32jl : 

R(21/31) = cot 2 fl 23 • (33) 

The results for the double ratio R(21/31) discussed above are illustrated in Fig.^ where 
the dependence of R(21/31) on the leptogenesis CPV parameter c for S13 = 0; 0.10; 0.20 and 
few characteristic values of the Majorana and Dirac CPV phases (a — (3m) = 0; 7r/2; ±7r and 
5 = 0; ±7r/2; ±7r are shown. The figure was obtained using the best fit values of the solar and 
atmospheric neutrino oscillation parameters #12, Am^, #23 and Am^. The lightest neutrino 
mass ttii was set to 0. The quantities |(Vjl^)f 2 H 13 | 2 were calculated using eqs. (JTHJ) and (|T7|) 
and not the approximate expressions given in eqs. (|20j) and (|21|). The leptogenesis CPV 
violating parameters a and b can contribute only to the higher order corrections 0(r 2 , s 2 3 ) 
in |(VjY" i/ )f 2 H 13 | 2 . These corrections can be relevant for the evaluation of R(21/31) in the 
case of cancellations between the terms in |(VjVj y )5 v 2 H 13 | 2 , which provide the leading order 
contributions. We have allowed a and b to vary in the same interval as the parameter c in 



11 



the calculations. The effects of the higher order corrections due to a and b is reflected in the 
widths of the lines in Fig. 1. 

We shall perform next similar analysis for the double ratio R(21/32) = BR(/i — > e + 
7)/BR(r — > fi + 7). As can be easily verified using eqs. (|2U|) and (j22j) and the known values 
of the neutrino oscillation parameters, for |c| < 0.3 we always have 

R(21/32) < 1 . (34) 

Typically the stronger inequality R(21/32) <C 1 holds 8 (see further). 

It is not difficult to convince oneself also that the term oc A 31 s 23 c 23 dominates in 
\(Y}Y U ) 23 \ 2 . Indeed, we have (Am^/Am 2 ^ ^ 0.18, c 12 |cos0 23 | ^ 0.24, s 13 Si 2 sin 2fl 23 ^ 
0.12, and for |c| < 0.2 (0.3), the terms oc c in eq. (J22|) give contributions which do not 
exceed approximately 8% (18%). Keeping only the largest of these contributions we have: 
|(n t n)2 3 | 2 = M\v-J |v^4 S23 c 2 3 + 2zc (Am2 1 Am| 1 )3c 12 cos(a - (3 M )/2\ 2 - Thus, the 
branching ratio BR(r — > ji + 7) exhibits very weak dependence on c and (a — (3m)- Up to 
the indicated corrections which for |c| < 0.3 can increase 1(1^1^)231 2 by not more than 18%, 
we have: 



M 



\{YlY v )^ = ^ (35) 



v 4 



Thus, for |c| ^ 0.3 in the case under discussion, BR(r — > fi + 7) depends essentially only on 
the atmospheric neutrino oscillation parameters Am^ and (9 23 (and not on the Dirac and 
Majorana CPV phases, leptogenesis CPV parameters or solar neutrino oscillation parameters 
AiUjj and 6* 12 ) and has a relatively simple form: BR(r — > /i+7) = Fx (Am 31 /(4t> 2 )) sin 2 2# 23 , 
where the factor F oc M|,/t> 2 contains all the dependence on M R , tan/3 and the SUSY 
breaking parameters (see eq. (JHJ)). The double ratio R(21/32), however, depends in the 
case under discussion both on c e % 2^ and Si 3 e _4<5 : 



\ r p nh -u a n NH r 
R(21/32) = |C23P A 2 + T? ' , (36) 



Am 31 s 23 c 23 



where P NH and Q NH are given by eqs. (}2Tj) and 

For si 3 = 0.2 and |c| ^ 0.25, we have \/ Am 31 si 3 = 2.3a/ Am^cuSu and a/ Am 31 si 3 ^ 
1.7 (2c(Am3 1 Am2 1 ) 3 Si 2 ). The double ratio R(21/32) exhibits noticeable dependence on the 
CPV phases (a — (3m) and 5. For (a — /3jv/) = and 5 = 0, the term oc c in Q NR gives a 
subdominant contribution and R(21/32) is practically independent of c. If 5 = tt, however, 
the term oc a/ Am^s^ in Q NH can be compensated partially by P NH and for sufficiently 
large values of |c| the term oc c in Q NH can be non-negligible. For |c| ^ 0.1 in this case we 
can have R(21/32) ~ few x 10~ 2 , while if (a - (3 M ) = tt and |c| ^ 0.2, R(21/32) can be as 
small as R(21 /32) ~ few x 10~ 3 . 

In the case of si 3 ~ \J Am^c^s^/ a/ Am^ ~ 0.07, partial compensation between the 
three terms in the numerator of the double ratio R(21/32) can take place. The double ratio 
R(21/32) can be particularly strongly suppressed for 5 = n, when values of R(21/32) ~ 

8 This is in contrast to the case of normal hierarchical heavy Majorana neutrino mass spectrum, in which 
one typically has R(21/32) - 1 [30]. 
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(10 3 — 10 4 ) for |c| ~ 0.05 are possible. Similar mutual compensations between the terms 
in the numerator of R(21/32) can be realised if S13 <C a/ Am^c^s^/ a/ Am 2 ^ an d \c\ ~ 
(0.15 - 0.20). One can have R(21/32) ~ (10^ 3 - 10~ 4 ) in this case as well. For sufficiently 
small S13 the dependence on the phase 5 is obviously insignificant and we have: 



R(21/32) 



\<m> 



NH 



a/A 



m 31 S 23 C 23 



Am^ 4 . . „ .- «-0m 



Am 2 31 



C23C12 + i 2c s 23 e* 



(37) 



If in addition 2\c\ <C (Am^/Amjj)* = 0.42, R(21/32) is also practically independent of c 
and (a — (3m)- It is determined completely by the solar and atmospheric neutrino oscillation 
parameters: 

R(21/32) = |<m>| NH 4 = ~ 1.3 x 10~ 2 . (38) 

31 23 31 23 

The specific features of the double ratio R(21/32) discussed above are evident in 
Fig. 2, where the dependence of R(21/32) on the leptogenesis CPV parameter c, |c| < 0.25, 
for three values of S13 = 0; 0.1; 0.2, and several characteristic values of the Major ana and 
Dirac CPV phases (a — (3m) and 5 is shown. Figure 2 was obtained using the same method 
and the same best fit values of the oscillation parameters Am^, sin 2 9\ 2 , Am^ and sin 2 2# 2 3, 
as Fig. 1. 



4.2 Inverted Hierarchical Spectrum 

If neutrino mass spectrum is inverted hierarchical (IH) one has 1713 <C mi^, and we 
shall assume that the terms oc ^fml in eqs. (j2~Uj) - (j2~2"j) can be neglected. For <m>in we 
have <m>iH = \l\Km^\{c\ 2 + s\ 2 e ta ) (see eq. ©). Now i 7^ j, depend on the 

Majorana phase a, on the leptogenesis CPV parameter a and, if S13 has a value close to the 
current upper limit - on the Dirac phase 5: 

Mr 



IH I 

12 I 



I D IH , ^vIHI 
C 23 P + S 2 3 Q 



\(yX)™\ 



M, 



R 



-s 23 P m + c 23 Q 



IHi 

23 I 



V 

u 

M R 



IHI 



I. . 



Am 2 



31 1 C 23 S 23 



a 



•1 + 4ac 12 Si2 sin — 



where 



jIH 



Am 2 x 



2 y/\Am 



C12S12 + i 2a <m> IH e" 



Q 



IH 



Am 2 



31 1 S 13 e 



a: 



1 + 4aci 2 Si 2 sin — 



For S13 satisfying 
sin 9 
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1 + 2|a| sin26» 12 ) < min 2|a| cos 29 



Ami 
7!2 '4|A^ Sin ^ 12 



(39) 
(40) 
(41) 

(42) 
(43) 

(44) 
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the dependence of lO^jKOif 13 1 on the Dirac phase 5 would be insignificant. The terms oc Q m 
in eqs. ff3~9*j) and (j4*Uj) are negligible, and the ratio of BR(/i — > e + 7) and BR(r — > e + 7) is 
given by 

R(21/31) = cot 2 fl 23 , (45) 

independently of the values of the Majorana CPV phase a, leptogenesis CPV parameter a, 
etc. (Fig. 3). If in addition \a\ < (Am| x /(8| Am 2 n \)) sin 26 12 ^ 3.6 x 10~ 3 , BR(/x -► e + 7) 
and BR(r — > e + 7) also will not depend on a and a: 

1(^^)11(13) I 2 = Ci2(i3) ( M R Am| 1 K)(Am 2 1 /(16|AmI 1 |)) sin 2 20 12 , where C 12(13) = c 23 (*«). 
In the case of \a\ cos26 , i 2 3> (Am| 1 /(8|Am| 1 |)) sin2#i 2 = 4 x 10~ 3 , however, we have: 

I(^)i!(i3)l = M \<m>m\ ^ C 12(13) . (46) 

u 

Thus, both BR(/i — > e + 7) and BR(r — > e + 7) are proportional to |a| 2 |<m> IH | 2 . 

We get R(21/31) ~ 1 also when s 13 > (Am 2 1 /(4| Am^D) sin 2# 12 = 8 x 10~ 3 , provided 
a = and 5 = 0; vr (Fig. 3). In this case |<m> IH | 2 = |Am^|, 1(1^1^)^1 = (4a 2 c 2 3 + 
S 2 3S 2 3 )|Am 2 1 |M 2 K, = (4a 2 , 2 3 + S 2 3 c 2 3 )|Am 2 1 |M 2 K and " 



R(21/31) = A / 2 C f + Si ff . (47) 
4a s 23 + s 13 c 23 



If, however, a is significantly different from zero, say a = ±7r/2; ± tt, and \a\ is sufficiently 
large, being comparable in magnitude to S13, the terms oc P IH and oc Q IH in |(^J^)if | or 
KyjY^)™! can partially compensate each other and we can have R(21/31) ~ (10~ 3 — 10~ 2 ) or 
R(21/31) ~ (10 2 — 10 3 ) (Fig. 3). For given \a\ and S13, the degree of compensation depends 
on the values of a and 5 and on the sgn(a). It is maximal in | (l^j"y^)™ ( - 13 - ) | , for, e.g., a = tt 
and 5 = 0, or a = —tt and 5 = tt, and a > (a < 0) (Fig. 3). 

The ratio of BR(/i — > e + 7) and BR(r — > \i + 7) depends both on a and a. For \a\ < 0.3 
we have R(21/32) ^ 1; if \a\ < 0.1, the stronger inequality R(21/32) <C 1 typically holds. 
For \a\ <C (Am2 1 /(8|Am 31 |)) sin2# 12 = 4 x 10~ 3 and negligibly small S13, for instance, one 
finds [30] R(21/32) = 10~ 4 . If, however, the term oc a dominates in \P |, i.e., if \a\ cos2# 12 ^> 
4 x 10~ 3 , we get (for S13 ~ 0) BR(« — > e + 7) oc \a\ 2 |<m>iH| 2 and correspondingly, 

R(21/32) = 4 |a| 2 S2 3 2 r m -1 + 2a V (1 - r ffl )^ * , (48) 
where 77 = sgn(sin 2#i2 sin |) and 

rm = = 1 - sin 2 29 12 sin 2 % . (49) 

|Am 31 | 2 

Now R(21/32) can be considerably larger: for a varying between and tt and \a\ having a 
value, e.g., in the interval (0.04 - 0.10), the ratio of interest satisfies 1.9 x 10~ 3 ^ R(21/32) ^ 
8.0 x 10" 2 , the maximal value corresponding to \a\ = 0.1 and a = 0. 

The predictions for the double ratios R(21/31) and R(21/32), corresponding to IH light 
neutrino mass spectrum are illustrated in Figs. 3 and 4, respectively. As in the case of Figs. 
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1 and 2, the quantities |(yjY^)™| 2 , i ^ j, have been calculated using eqs. (JT5]) and (|T7|) 
rather than the approximate expressions given in eqs. ()20|) and (}2~T]) . The lightest neutrino 
mass m3 set to 0. The leptogenesis CPV parameters b and c, which can contribute only to 
the higher order corrections in |(yjy^)iH| 2 Q f interest, were varied in the same interval as the 
parameter a in the calculations. The effects of the higher order corrections due to b and c is 
reflected in the widths of the lines in Figs. 3 and 4. 



4.3 Quasi-Degenerate Neutrinos 



For QD light neutrino mass spectrum, m.1^,3 — m ^ 0.1 eV, one has <m>Q£, 
m ( c i2 + s i2 e * a )' an d y^mj = m in eqs. (J2U|l - (j22)) - Barring "accidental" cancellations, we 
always have |(FjYy)i 2 | ~ |(Y^Y„)i3|, and correspondingly BR(jU — > e+7) ~ BR(r — - > e+7), in 
this case (Fig. 5). The expressions for | (YjY u ) 12(13) | of interest simplify if max(|a|, |6|, |c|) 
max(Am2i/(4m 2 ), | Am3 1 |si 3 /(4m 2 )) and S13 ^ 0.1: 

\(YX)??\ = 2 ^ |c 23 P QD + s 23 Q^\ , (50) 

v 

u 

\(yX)?s\ = ^ |- S 23P QD + c 23 g QD | , (51) 

u 

where 

P QD = a<m> QD e- i f , (52) 
gQ° = m (bcx2 + cs 12 e^) + i aSl3 sin 26> 12 e~ iS sin ^ . (53) 

The condition specified above is compatible with the leptogenesis constraints on the product 
\abc\ [26]. For a 9* and (3 M = ±vr we get: 

2 IpQDI 2 , 2 InQDI 2 

R(21/31) = 2 2 3 pqd 2 f XJ2 • (54) 
s 2 3 |PQ D | + c 2 3 |QQ D | 

Obviously, in this case either R(21/31) = 1 independently of the value of #23, or R(21/31) = 
tan 2 #23 or cot 2 #23. 

Under the condition leading to eqs. (joHj) - (|S3*j). the quantity |(FjF„)23|, eq. (|2"2"j). cannot 
be simplified. The term oc A 3i in eq. (|22j) will be the dominant one if 
max(|asin(a/2)|, 2|a|si 3 , l a | 2 ), \b\, l c l ^ |Am3 1 |/(4m 2 ). Given the leptogenesis constraint on 
\abc\, this is realised, e.g., for m = 0.1 eV if \a\ ~ |6| ~ |c| = 10~ 2 , or if sin(a/2) = 0, S13 = 
and \a\ 3> |6|, |c| but \a\ 2 <C | Am 2 1 |/(4m 2 ). In both cases we have 

2 

< 1 . (55) 



R(21/32)- 16m4 kl 



(Am 



31/ 



<m> QD 6ci2 + csi 2 e*2 _ i ^ MZ 
— I e 2 



m s 23 a c 23 



For m = 0.10 eV and \a\ = \b\ = \c\ = 10 2 , the ratio R(21/32) given by eq. (|B3j) depends 
on a, Pm, sgn(6/a) and sgn(c/a) and satisfies 2 x 10 -4 ^ R(21/32) ^3x 10 _1 . If, however, 
a S 0, si 3 = and |a| = 0.2 with |a6c| = 10" 5 , the "corrections" oc |a| 2 in eq. (J22J) will 
be non-negligible since |a| 2 ~ |Am 2 1 |/(4m 2 ). In this case we can have even R(21/32) = 200 
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as a consequence of rather strong partial cancellation between the different terms in the 
expression for | (1^1^)23! (Fig- 6). 

The term oc A31 in eq. (j22j) can be neglected if, e.g., at least one of the CPV param- 
eters \a sin(a/2)|, \b sm((3 M /2) cos 26*23! (1^ cos (/^m/2)| 2 ) and \c sin((a — (3m)/2) cos20 2 3| 
(|c cos((a — /?m)/2)| 2 ) is much bigger than jAm^l/ '(4m 2 ) (|Am 2 1 | 2 /(4m 2 ) 2 ). In this case 
eqs. (J5UJ) - (JS2I) are also valid. We get particularly simple expressions for |(Yj"Y,,)^ D |, i 7^ j, 
if the terms oc a in eqs. pUj) - (f2~2~j) dominate: 

\(YX)??\ = 2 |a| ^ |<m>| QD c 23 , (56) 
K y jK)? 3 D | = |(^)? 2 D | tan ^3 , (57) 
|(yJn) 2 Q 3 D | = 2 |a| ^ v /m2-|<m>| 2 1D c 23 s 2 3 . (58) 

Equations (|5^j) - (J57)l are valid provided \a\ 3> max(|6|, |c|, | Am 2 1 |/(4m 2 )), while eq. ([HSjl 
holds if |a sin(a/2)| 3> max(|6|, |c|, | Am 2 1 |/(4m 2 )). For \a\ < 1 and, e.g., m = 0.1 eV, the 
latter condition requires | sin(a/2)| = 1. In these cases both BR(/i — > e + 7) ~ |a| 2 |<m>|Q D 



and BR(r — > e + 7) ~ |a| 2 |<m>|Q D , while BR(r — > /i + 7) ~ |a| 2 (m 2 — |<m>|Q D 
a| 2 m 2 sin 2 2#i2 sin 2 (a/2). For the ratio of the first two we get 

R(21/31) 2* cot 2 #23 , (59) 



which should be compared with eqs. (j2HJ) and (j4^|) . The ratio of BR(/i — > e + 7) and 
BR(r — > + 7) is independent of the leptogenesis CPV parameter a. Given 6 I 12 and 6*23, it 
is determined by the Majorana phase a: 



R(21/32) ~ , - ? Q 2 , , = 2 12 2 ■ (60) 





<m> 


qd _ 1 -sin 2 2fl 12 sin 2 (a/2) 


(m 2 — 


<m>\ 


QD ) ^3 s 2 3 sin 2 2# 12 sin 2 («/2) 



We get similar results if the terms oc b (oc c) dominate in eqs. (}2T)j) - (}2*2*j) . which in the 
case of eq. (|2^j) would be possible only if the Majorana phase (3m (phase difference a — (3m) 
deviates significantly from tt. Now K^jY"^)^! depends on (3m (a — (3m)- If, e.g, the terms 
oc c dominate we get: R(21/32) = s^tan 2 ^ (1 - sin 2 2# 23 sin 2 ((a - p M )/2)y 1 . 
Our results for the double ratios R(21/31) and R(21/32) are illustrated in Fig. 5. 

5 Conclusions 

Working in the framework of the class of SUSY theories with see-saw mechanism and 
soft SUSY breaking with flavour-universal boundary conditions at a scale M x > M R , we 
have analysed the dependence of the rates of lepton flavour violating (LFV) decays /i — > e+7, 
t — > e + 7 r — > fi + 7 (li — > l 3 • + 7) and of their ratios, on the Majorana and Dirac CP- 
violation (CPV) phases in the PMNS matrix ?7pmns, ®, (3m and 5, and on the leptogenesis 
CP-violating (CPV) parameters. The case of quasi-degenerate in mass heavy RH neutrinos 
was investigated, Mi = M 2 = M 3 = Mr, assuming that splitting between the masses of the 
heavy neutrinos is sufficiently small, so that it has practically no effect on the predictions for 
the U — > lj + 7 decay rates. Results for the normal hierarchical (NH), inverted hierarchical 
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(IH) and quasi-degenerate (QD) light neutrino mass spectra have been derived. The analysis 
was performed under the condition of negligible renormalization group (RG) effects for the 
light neutrino masses rrij and the mixing angles and CPV phases in ?7pmns 9 - in the wide 
region of validity of eqs. and (|13|) in the relevant SUSY parameter space, the ratios of 
rates of the decays ji — > e + 7 and r — > e + 7, and \i — > e+7 and r — > /i+7, are independent of 
the SUSY parameters - they are determined by the neutrino masses (rrij) and mixing angles, 
Majorana and Dirac CPV phases and by the leptogenesis CPV parameter (s). For the matrix 
of neutrino Yukawa couplings, Y u - a basic quantity in the analysis performed, we have used 
the orthogonal parametrisation [23]. The latter proved to be the most convenient for the 
purposes of our study [26]. In this parametrisation Y u is expressed in terms of the light and 
heavy Majorana neutrino masses, £7pmns, and an orthogonal matrix R [23]. Leptogenesis can 
take place only if R ^ 1 (see, e.g., [35,38]). In the case of quasi-degenerate in mass heavy 
Majorana neutrinos considered in the present article, the rates of the LFV decays k — > lj +7 
of interest do not depend on the matrix R ^ 1 if R is real. For complex matrix R, only 
the three leptogenesis CPV real dimensionless parameters of R, a, b and c (eqs. (|17|) and 
(HHJ)), enter into the expressions for the k — > lj + 7 decay branching ratios of interest [26], 
BR(/j — > lj + 7). In our analysis we have assumed that |a|, |c| < 1, as is suggested by the 
leptogenesis constraint derived for QD light neutrinos and negligible favour effects [26]. In 
various estimates we have considered values of \a\, \b\, \c\ < 0.3. The case of real matrix R 
corresponds effectively to a,b,c = 0. 

We have found that for NH (IH) spectrum and negligible lightest neutrino mass mi 
(7713), the branching ratios BR(/j — > lj + 7) depend, in general, on one Majorana and the 
Dirac CPV phases, a — (3m (at) and 5, one leptogenesis CPV parameter, c (a), on the 
CHOOZ angle #13 and on the mixing angles and mass squared differences associated with 
solar and atmospheric neutrino oscillations, # 12 , Am^, and 823, Am^. The double ratios 
R(21/31) oc BR(/i -> e + 7)/BR(r -> e + 7) and R(21/32) oc BR(/x -> e + 7)/BR(r -> /1 + 7) 
(see eq. (J2U0 are determined by these parameters. The same Majorana phase a — (3m 
(a) enters also into the NH (IH) expression for the effective Majorana mass in neutrinoless 
double beta ((/3/3)o u -) decay, <m> (eqs. (JSJ) and ©)• For the QD spectrum, BR(/j — > lj+j) 
depend, in general, on the absolute neutrino mass m, the three leptogenesis CPV parameters, 
a, b, c and on the two Majorana phases a and (3m- For the IH and QD spectra, the phase 
a enters into the expressions for BR(/i(r) — > e + 7), in particular, through the effective 
Majorana mass <m> (see eqs. (jj2j) and (|52|l). Our results for the double ratios show that 
we can have R(21/31) ~ 1 or R(21/31) < 1, or else R(21/31) > 1 in the cases of NH and IH 
spectra, while for the QD spectrum typically R(21/31) ~ 1. In contrast, for the NH and IH 
spectra one always gets R(21/32) < 1; in most of the relevant parameter space R(21/32) 1 
holds. For the QD spectrum, however, R(21/32) ^ 1 is also possible. 

More specifically, we find that for the NH (IH) spectrum, BR(ju(r) — > e + 7) exhibit 
significant dependence on the leptogenesis CPV parameter c (a) and on the Majorana CPV 
phase a — (3m («) for |c| ?Z 0.02 (|a| ^ 0.02) and for any S13 ^ 0.1 (si 3 ^ 0.2). In 
certain cases the dependence of BR(/x(r) — » e + 7) on the phase a — (3m (of) and/or the 
parameter c (a) is dramatic. More generally, the dependence of BR(/i(r) — > e + 7) on 

9 It is well-known that in the class of SUSY theories considered, this condition is satisfied in the cases of 
NH and IH light neutrino mass spectra; it is fulfilled for the QD spectrum provided the SUSY parameter 
tan/? is relatively small, tan/3 < 10. 
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the Majorana phase can be noticeable only if the corresponding leptogenesis parameter is 
sufficiently large: for \c\ <C max((Am| 1 /Am2 1 ) 3 Si3, 0.5(Am| 1 /Am| 1 )4) in the NH case, 
and \a\ <C max(si 3 / (2 cos 2#i2), Am^ tan 2^12 / (8 1 Am| x | ) in the IH one, both c (a) and the 
Majorana phase have practically no effect on BR(/x(r) — > e + 7). Similarly, the CHOOZ 
angle #13 and the Dirac phase 5 can be relevant in the evaluation of BR(/i(r) — > e + 7) 
in the cases of NH and IH spectra only if s 13 is large enough, i.e., if respectively, s 13 <; 
^Am 2 21 sin20i 2 /(2v'Am§ 1 |) = 0.07, and s 13 > Am| 1 sin20 a2 /(2|Am§ 1 |) = 8 x 10" 3 . In the 
case of NH (IH) spectrum, BR(r — > \i + 7) is practically independent of si 3 ^ 0.2; the 
dependence of BR(r — > \i + 7) on the leptogenesis parameter c (a) and the Majorana phase 
a — Pm (a) is relatively weak for |c| ^ 0.3 (\a\ ^ 0.1). For this wide range of values of |c| 
(|o|) we have BR(r -> // + 7) = F x (| Am^|/(4^)) sin 2 2# 23 , where the factor F oc M 2 R /v 2 
contains all the dependence on Mr, tan/3 and the SUSY breaking parameters (see eq. (fTTj)). 

The double ratios R(21/31) and R(21/32) (Figs. 1-4) can exhibit in the cases of NH and 
IH spectra strong dependence on the Dirac and/or Majorana phases if S13 ~ 0.1 — 0.2 and/or 
if the relevant leptogenesis parameter exceeds approximately 10~ 2 . Under the indicated 
conditions values of R(21/31) ~ (10~ 3 - 10~ 2 ) < 1 or R(21/31) ~ (10 3 - 10 2 ) > 1, are 
possible. For, e.g., si 3 ~ 0.1, the sign of the inequality is determined by the sign of the 
leptogenesis parameter, the value of the Majorana phase and/or the value of the Dirac phase 
(Figs. 1 and 3). If for the NH (IH) spectrum, a - f3 M = (a = vr) and 5 = ±tt/2, ±3tt/2, 
R(21/31) takes one of the following three values R(21/31) = 1; tan 2 6> 23 ; cot 2 6» 23 . For \a\ > 
Am| 1 tan26'i 2 /(8|Am^ 1 |) in the IH case, we find BR(/x(r) — > e+7) = F IH x \a\ 2 \<m>m\ 2 /v 2 ., 
and thus R(21/31) = 1, where <m>in is the effective Majorana mass in (/3/?) 0l/ -decay and 
the factor F IH oc M^/v 2 includes the dependence on M R and on the SUSY parameters. 
For sufficiently small s i3 and |c| (s 13 0.07, |c| 0.2)) in the case of NH spectrum, we 
get: R(21/32) ^ Am 2 1 /(Am 2 1 s^ 3 ) c\ 2 s\ 2 = 10" 2 . Smaller values of R(21/32) are possible, 
e.g., for S13 = (0.1 — 0.2), if |c| ~ 0.05 and if for given sgn(c), the Majorana and Dirac 
phases (a — (3m) and 5 have specific values (Fig. 2). For the IH spectrum we typically have 
R(21/32) < 1 for \a\ < 0.1. If 2\a\ cos2#i 2 , sin^i 3 < 0.5Am 2 21 c 12 s 12 /\ Am^l, R(21/32) is 
completely determined by the solar and atmospheric neutrino oscillation parameters Am 21 , 
012, Am 2 x and 8 23 , and R(21/32) ^ 10~ 4 . 

In the case of QD light neutrino mass spectrum, the leptogenesis constraint implies [26] 
10~ 6 ^ \abc\ < 10~ 4 . The expressions for BR(Zj — > lj + 7) and for the double ratios R(21/31) 
and R(21/32) simplify considerably if the terms including one given leptogenesis parameter 
dominate. We get, e.g., R(21/31) = tan 2 #23 and R(21/32) = 1 if the terms oc b (oc c) 
are the dominant one. This requires relatively large values of \a\ or |6| or |c|. If, however, 
\a\ ~ |6| ~ |c| ~ 10^ 2 , R(21/32) lies in the interval ~ (10~ 4 - 10" 1 ). 
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Note Added. During the completion of the present study we became aware [58] that 
an analysis along seemingly similar lines is being performed by R. Riickl et al. 
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Figure 1: The double ratio R(21/31) in the case of NH light neutrino mass spectrum, as a 
function of the leptogenesis CPV parameter c, for s 13 = 0.2; 0.1; and several characteristic 
values of the Dirac and Majorana CPV phases 5 and a — (3m- The figure was obtained using 
the best fit values of the solar and atmospheric neutrino oscillation parameters Am^, sin 2 #12, 
Amjj and sin 2 2#23- The lightest neutrino mass m\ was set to 0. The effects of the higher 
order corrections in leptogenesis CPV parameters is reflected in the width of the lines (see 
text for further details). 
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Figure 2: The same as in Fig. [TJ but for the double ratio R(21/32). 
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(e) (f) 

Figure 3: The double ratio R(21/31) in the case of IH light neutrino mass spectrum as a 
function of the leptogenesis CPV parameter a, for several characteristic values of the CHOOZ 
angle #13 and Majorana and Dirac CPV phases a and S. The results shown correspond to 
the lightest neutrino mass m 3 = 0. The effects of the higher order corrections due to the 
leptogenesis CPV parameters b and c is reflected in the widths of the lines (see text for 
further details). 
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Figure 5: The double ratios R(21/31) and R(21/32) as a function of the leptogenesis CPV 
parameter a in the case of QD light neutrino mass spectrum for S13 = 0;0.1 and several 
values of the Majorana phases a and (5m- The results shown are obtained for \abc\ = 10~ 5 
and \b\ = |c|. For s 13 = 0.1, values of the Dirac CPV phase < 5 < ti were considered. The 
lightest neutrino mass is set to m\ = 0.1 eV. 
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